Ultra-cold fermions in the flatland: evolution from BCS to Bose 
superfluidity in two-dimensions with spin-orbit and Zeeman fields 
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We discuss the evolution from BCS to Bose superfluidity for ultracold fermions in two-dimensions 
and in the presence of simultaneous spin-orbit and Zeeman fields. We analyze several thermodynamic 
properties to characterize different superfluid phases including pressure, compressibility, induced 
polarization, and spin susceptibility. Furthermore, we compute the momentum distribution and 
construct topological invariants for each of the superfluid phases. 
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PACS numbers: 03.75.Ss, 67.85.Lm, 67.85.-d 

Ultra-cold Fermi atoms are one of the most interesting 
physical systems of the last decade, as they have served as 
quantum simulators of crossover phenomena and phase 
transitions encountered in several areas of physics. Due 
to their tunable interactions, atoms like 6 Li and 40 K have 
been used to study the crossover from BCS to BEC super- 
fluidity, to simulate superfluidity in neutron stars, and to 
investigate unitary interactions which are of great inter- 
est in nuclear physics. Furthermore, the ability to con- 
trol the internal spin state of the atoms by using radio- 
frequencies (RF) enabled the studies of quantum and 
classical phase transitions as a function of interactions 
and population imbalance. These tools have permitted 
the study of crossover phenomena and phase transitions, 
and have validated the symmetry based classification of 
phase transitions put forth by Landau over the thermo- 
dynamic classification proposed earlier by Ehrenfest. 

Very recently three new tools have been developed for 
ultra-cold fermions. The first tool is a method that allows 
for the extraction of thermodynamic properties - such as 
pressure, entropy, compressibility and spin-susceptibility 
- for uniform systems by using local density profiles [lj, [2J . 
The second tool provides the ability to trap two dimen- 
sional fermions in the quantum degenerate regime 0, Q] ■ 
The third tool is the creation of artificial spin-orbit and 
Zeeman fields (H, which allows for the visitation of pre- 
viously inaccessible physical regimes, such as the limit 
of strong spin-orbit coupling (SOC). Thus far, the SOC 
created is of the equal Rashba-Dresselhaus (ERD) type 
instead of Rashba-only Q or Dresselhaus-only . The 
emergence of the third tool in the context of ultra-cold 
bosons [j| has generated substantial theoretical interest 
about the possibility of applying_the techni que to thrce- 
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dimensional ultra-cold fermions 
groups studied the Rashba SOC 
occurs naturally in condensed matter physics, while our 
group explored the ERD case 14, T^. In the context of 
ultra-cold fermions, the combination of the first two new 
tools has been achieved experimentally very recently fli| , 
and there are recent reports that ERD SOC has been cre- 
ated in fermionic systems of 40 K [TtJ and 6 Li 

In this manuscript, we discuss the combination of these 



three new techniques to describe the evolution from BCS 
to Bose superfluidity in two-dimensions, which in the con- 
text of ultra-cold fermions has been discussed for popula- 
tion balanced s-wave [3] , p-wave [2(3] , and unbalanced s- 
wave [2l| systems without SOC. Preliminary accounts of 



the effects of Rashba SOC have recently appeared in the 
cold atom literature [22|,[23j]. Here, we study the simulta- 
neous effects of ERD spin-orbit and Zeeman fields, which 
produce a very rich phase diagram consisting of different 
topological superfluid. To describe these novel phases, 
we discuss their spectroscopic and thermodynamic prop- 
erties. 

Hamiltonian: We start with the total Hamiltonian 
H = H sp + Hi n t, where the single-particle part is (?i = 1) 



H sp — 



J2 V'l(k) \l<5 ss , -Ai(kK, 



Mk). (1) 



k,s,s' 



Here, K = k 2 /(2m) — [i is the kinetic energy in reference 
to the chemical potential fj, and hi(k) is the spin-orbit 
plus Zeeman field along the /-direction (s,s' ='\\-l- i = 
x,y,z). In the present case, we consider only the ERD 
SOC field h ERD = (0,1^,0) [1,0 [H and Zeeman field 
hzEE = (0, h y ,h z ). The current Hamiltonian matrix that 
can be created in the laboratory using Raman beams 
is of the form -ffzso = —h z a z — h y a y — ZiERDCy, where 
h z = (0/2) with £1 being the Raman intensity, h y = (6/2) 
with 6 being the detuning, and /ierd = vk x with v being 
a measure of the SOC strength. 
The interaction Hamiltonian is 



H int = VW&i(k)6 q (k'), 



(2) 



k,k',q 



where 6 q (k) = ^{(k + q/2)^|(-k + q/2) is the creation 
operator for a pair of fermions with center of mass mo- 
mentum q and relative momentum 2k. The interaction 
term Vkk' is the Fourier transform of the contact inter- 
action V(r,r') = -g6(r - r') (g > 0). 

The single-particle Hamiltonian H sp can be diagonal- 
ized by a momentum dependent SU(2) transformation 
into the hclicity basis $t(k) = M (k)^^ (k) , where the 
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helicity creation operators $t(k) = f<^|(k), </>|(k)J are 
expressed as linear combinations of the standard cre- 
ation operators >]>t(k) = ^V{(k), (k)) ■ The eigenval- 
ues of the Hamiltonian matrix H sp (k) are ^^(k) = 
A'(k) =F |h ff (k)|, corresponding to the helicity spin ft (4) 
being aligned (anti-aligned) with respect to the effective 
magnetic field h c ff (k) = (0, /iERD(k), h z ). Here, the mag- 
nitude of the effective field is |h c ff(k)| = \Jh 2 + v 2 k 2 and 
A(k) = e k - (i = k 2 /(2m) - /a. 

In the helicity basis, the interaction Hamiltonian 

transforms as H int = -gJ2q a pj8 B if}( c l) B i^ c 0^ where 
a, /3, 7, (5 =ft, Pairing is now described by the operator 
B Ui<i) = £k A «£(ki,k 2 )$t (ki)$],(k 2 ) and its Hcrmi- 
tian conjugate, with momentum indices ki = k + q/2 
and k 2 = k + q/2. Pairing occurs between fermions 
of momenta ki and k 2 in two different helicity bands 
(inter-helicity pairing), as well as within the same helic- 
ity band (intra- helicity pairing). The tensor A Q( g(ki,k2) 
contains the matrix elements of the momentum depen- 
dent SU(2) rotation into the helicity basis, and reveals 
that the center of mass momentum ki + k 2 = q and the 
relative momentum ki — k 2 = 2k arc coupled and no 
longer independent. 

For pairing at q = 0, the order parameter for 
superfluidity is the tensor A Q( g(k) = AoA Q/ 3(k, — k), 
where A = — gJ2-ys( B -ys(0)), leading to components: 
A-fl-fi-(k) = iAx(k) for helicity projection A = +1; 
A^jj.(k) = — A^(k) = — Ag(k) for helicity projection 
A = 0; and Ajj.^(k) = — iAy(k) for helicity projection 
A = -1. The amplitudes A T (k) = A /i£flD(k)/|h eff (k)| 
and Ag(k) = Ao/i z /|h e fj(k)| reflect the triplet and sin- 
glet components of the order parameter in the helicity 
basis. For the ERD case, A^(k) is odd and A$(k) is 
even under parity V in momentum space, while they are 
both invariant under time reversal T, and both reflect 
the broken U(I) symmetry of the superfluid phase. The 
triplet and singlet sectors in the helicity basis are not in- 
dependent as the interactions from which they originate 
occur only in the singlet s-wave of the original spin basis, 
and the relation |A T (k)| 2 + |A s (k)| 2 = |A | 2 holds. 

It is worth emphasizing that in the triplet sector 
A^(k) and A,p.(k) contain not only p-wave, but also 
/-wave and higher odd partial waves in two dimensions, 
as can be seen from a multipole expansion of | h c ff (k) | _ 1 = 

2/ 2 

[h 2 + ^erdC 4 )] f° r unite h z . Similarly in the singlet 
sector Af^k) and Ajj.^k) contain s-wave, d-wave and 
higher even partial waves, as long as the Zeeman field h z 
is non-zero. Higher angular momentum pairing occurs 
because the local (zero-ranged) interaction in the (t, I) 
spin basis is transformed into a finite-ranged anisotropic 
interaction in the helicity basis (ft,-IJ-)- 

The effective Hamiltonian in the helicity basis takes 



the matrix form 



ffeff(k) 











%(k) A M (k)^ 
A M (k) A M (k) 



A;.(k) At.(k) 





V A; 4 (k) A* 4 (k) o -^(k) ) 

The eigenvalues for the quasiparticle bands arc 



-£t(k) 




E p± (k) = J(E s (k) ± |h cff (k)|r + |A T (k)|= 



(3) 



(4) 



where As(k) = y |A'(k)| 2 + |As(k)| 2 is a characteris- 
tic energy for the singlet sector. The eigenvalues for 
quasihole bands are A^±(k) = — E p ±(k). The effective 
field | h c g (k) | emerges from the energy difference between 
the helicity energies |h c ff (k)| = [£y.(k) — £^(k)] /2, while 
A(k) = [^(k) +^(k)] /2 is the average energy of the 
helicity bands. 

Only A p _(k) has zeros or nodes, provided that the fol- 
lowing conditions are satisfied simultaneously: first, the 
effective magnetic field energy must equal to the excita- 
tion energy for the singlet sector |h c ff (k)| = As(k); sec- 
ond, the triplet component of the order parameter must 
vanish | Ay(k)| = 0. The nodal structure for A p _ (k) dom- 
inates the physics at low energies. For the ERD case, 
|/iERD(k)| = i>|fc;c|- Therefore, zeros of A p _(k) must oc- 
cur when k x = 0, which gives rise to up to four possible 
point nodes that are determined by the real solutions of 
k y to the equation [k 2 /(2m) — /i) 2 + |A | 2 = h 2 . 

To obtain the phase diagram, we calculate the ther- 
modynamic potential Q = — Tin [Tre _ff ° ff / T ] , and 
solve sclf-consistcntly the order parameter equation 
<9f2/<9|Ao| 2 = and the number equation TV = 
— (dft/dfi) T v . The contact interaction g/V is expressed 
in terms of the two-body binding energy Eb in the ab- 
sence of SOC via the relation V/g = J2k l/(2ek + Eb). 

In Fig. [I] we show the zero temperature phase dia- 
grams of Zeeman field h z versus the two-body binding 
energy Eb as well as induced population imbalance Pi n( j 
versus binding energy for both vanishing and large finite 
SOC (v/vp = and 0.8). We label the uniform superfluid 
phases with zero, one or two pairs of nodes as US-0, US-1, 
and US-2, respectively. The US-2/US-I phase boundary 
is determined by the condition [i = ^Jh 2 — | Ao| 2 , when 
\h z \ > |A |; the US-O/US-2 boundary is determined by 
the Clogston-like condition \h z \ = |Ao| when fi > 0, 
where the gapped US-0 phase disappears leading to the 
gapless US-2 phase; and the US-0/US-I phase boundary 
is determined by \i = — \A 2 — |Ao| 2 , when \h z \ > |Ao|. 
Furthermore, within the US-0 phase, a crossover line be- 
tween an indirectly gapped and a directly gapped US-0 
phase occurs at /i = 0. Non-uniform (NU) phases also 
emerge in regions where uniform phases are thermody- 
namically unstable for v/vp ^= 0, however NU phases 
only exist in the regime ofO < v/vp < {v/vp) c , where 
(v/vp) c « 0.6. Possible NU phases include phase separa- 
tion, modulated superfluids or supersolid. These possible 
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FIG. 1: Zero temperature (T = 0) phase diagrams of Zee- 
man field hz/eF versus two-body binding energy Eb/eF for 
equal Rashba-Dresselhaus (ERD) SOC (a) v/vf = and (c) 
v/vf — 0.8. The corresponding phase diagrams for the in- 
duced population imbalance versus Ef,/eF are shown in 
(b) v/vf = and (d) v/vf = 0.8. Here, cf = kp/(2m) and 
Vf = kp/rn are the Fermi energy and momentum respectively, 
and kF is related to the particle density n via k F = (3?r 2 n) 1/3 . 



phases are not illustrated in Fig. Q] since we are interested 
in comparisons between the zero SOC case v/vf = and 
the large SOC case v/vp > (v/vf) c - In the latter case, 
a tri-critical point occurs at the intersection of the US-0, 
US-1, and US-2 (where fj, = and \h z \ = |A |). 
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FIG. 2: The order parameter amplitude |Ao|/cf as a 
function of external Zeeman field h z /tF for (a) v/vf = 0, 
T/e F = 0, (b) v/vf = 0, T/e F = 0.03, (c) v/v F = 0.8, 
T/e F = 0, (d) v/vp = 0.8, T/e F = 0.03 at various binding 
energies Eb/eF = 0.1 (solid line), Eb/eF = 0.2 (dashed line), 
Eb/eF = 0.5 (dotted line), E b /eF = 1.0 (dot-dashed line). 

It is important to note from the phase diagrams that 
for finite SOC (v/vf > 0) it is always possible to form 



pairs in the lower helicity band (k) no matter how large 
the Zeeman field h z is. This is because an induced triplet 
component of the order parameter emerges and circum- 
vents the standard pair breaking Clogston limit for sin- 
glet pairing. Thus, the stable superfluid phase for large 
h z at T = is the US-1 phase as shown in Fig.[T]providcd 
that v/vp 7^ 0. This is clearly seen in Fig. [2] where the 
order parameter is shown for the cases of v/vf = and 
v/vf = 0.8. In the first case (v/vf = 0), the Clogston 
limit at T = occurs when h z /tF ~ |Ao|/cf as seen in 
Fig. where the order parameter jumps discontinously 
to zero at a critical value of (h z /ep)c- The discontinuity 
occurs also at finite temperatures (Fig.[2}D), rendering the 
transition from the superfluid to the normal phase discon- 
tinuous according to Landau's classification or first order 
according to Ehrenfest's. For v/vp = 0.8, the absence of 
the Clogston paramagnetic limit at zero temperature is 
clearly seen in Fig. [2}:, where the order parameter van- 
ishes only asymptotically as h z /ep — > oo because of the 
induced triplet component in the helicity basis. As the 
temperature is raised, thermal fluctuations break pairs, 
the order parameter vanishes, and the normal state is 
reached beyond a critical Zeeman field (h z /ep)c, which 
is an increasing function of the binding energy Eb /e f for 
fixed v/vf > {v/vp) c - This also occurs for Rashba-only 
couplings in two-dimensions [221 ]. and for ERD couplings 
in three-dimensions [lij . 

These superfluid phases can be classified according to 
the topological properties of their low-energy excitation 
spectrum £J p _(k). The nodes of i? p _(k) can be viewed 
as vortex singularities in momentum space. Consider the 
two-dimensional unit vector rii(k) = (m m ,m 



y ), where 



c(k) 



[£ s (k)-|h eff (k)|]/£ p _(k) and m y (k) 



Ar(k)/ E p -(k). The topological charge is the winding 
(Chern) number N w — (27r) _1 § d£z ■ m x dm/d£, where 
the line integral is over any line in the k x -k y plane enclos- 
ing a zero of E p _ (k). The topological charges take values 
±1. In the US-2 (US-1) phase the topological charges fol- 
low the sequence H h— (H — ) for the four (two) nodal 

points as we scan k y from negative to postive. 

The change in nodal structure is associated with bulk 
topological transitions of the Lifshitz class as noted for 
p-wave (2^, 24 1 and <i-wave [25|, [2(| superfluids. The 
loss of nodal regions of E p _ (k) correspond to annihila- 
tion of Dirac quasi-particles with opposite momenta and 
opposite topological charge (US-l/US-2 and US-l/d-US- 
0), or with the same momenta but opposite topological 
charge (US-2/i-US-0). For instance, starting from the 
US-1 (US-2) phase, Dirac quasi-particles annihilate at 
zero momentum producing bulk Majorana quasi-particles 
at the US-1 /d-US-0 (US-2/US-1) phase boundary, which 
then become massive Dirac fcrmions near zero momen- 
tum with positive mass (negative mass) in the d-US-0 
(US-1) phase. However, zero-momentum bulk Majorana 
quasi-particles do not occur at the US-2/i-US-0 phase 
boundary, as the four node superfluid US-2 leads be- 
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FIG. 3: (color online) The momentum distributions 
n a {k x ,k y ) for ERD SOC v/v F = 0.8 and E b /e F = 0.1 at 
T = 0, where s ="\ (\) for upper (lower) panels. (a)(d) i-US-0 
phase with h z /tF = 0.2; (b)(e) US-2 phase with h z /eF = 0.4; 
(c)(f) US-1 phase with h z /eF = 1-0. The color coding varies 
continuously from purple (n s = 0) to red (n s = 1). 



comes a doubly-degenerate two-node superfluid at finite 
momentum at the phase boundary and then becomes an 
indirect gap superfluid. 

The emergence of Dirac fcrmions produces dramatic 
changes in momentum distributions [25| 



l-5>^(k)] 



dfi s 



(5) 



which arc illustrated in Fig.[3]for parameters v/vp = 0.8 
and Eb/tp = 0.1 at T = 0. Different topological phases 
are shown for different values of Zecman field h z . Here, 
n F [Ej(]c)} = l/(exp[.Ej'(k)/T] + 1) is the Fermi function. 
For the i-US-0 phase the momentum distributions are 
smooth functions of k = (k Xl k y ) until the i-US-0/US- 
2 phase boundary is reached where doubly degenerate 
nodes emerge at k = (0, ± v / 2m/i). Beyond this point 
discontinuities develop in n^ s , which is due to the fi- 
nite Zeeman field and can be best seen in its f com- 
ponent. Within the US-2 phase there are four discon- 
tinuities at the zeros of ^-(k), which are located at 



0,±W2m 



Around these ze- 



ros the energy has cusps and the dispersion is linear. 
As h z is increased further, the US-2/US-1 boundary is 
crossed and the inner Dirac fermions annihilate, leaving 
discontinuities in the momentum distribution occur only 



at two nodes k 



0,±W2m 



The 



final US-1 phase is highly polarized as can be seen from 
the great contrast between Fig. [3fc and |3f . 

We show in Fig. |4] several thermodynamic properties 
as a function of Zeeman field h z jtF for different binding 
energies at T = 0. In Fig. IJJ,, we show the pressure P = 




FIG. 4: Various thermodynamic quantities versus Zeeman 
field h z /eF for different values of binding energies Eb/tF = 0.1 
(solid line), E b /e F = 0.2 (dashed line), E b /e F = 0.5 (dotted 
line), and Eb/e F = 1.0 (dot-dashed line), (a) The pressure 
P; (b) The isothermal compressibility N 2 kt\ (c) The induced 
polarization Pi n d; (d) The spin susceptibility \ zz . 



— fl/V. In Fig. [4Jd, we show the isothermal compressibil- 
ity kt = —V~ 1 (dV/dP)T, which can also be written as 
k t = N^ 2 (dN/d^) T - In Fig. Hh and SH, we show the 
induced polarization Pi nc j = (Nf — N±) / (N^ + N±), and 
the spin susceptibility \zz = (dPmd/dh z ) T v . Accord- 
ing to Fig. QJ>, if we fix the value of Eb and increase the 
Zeeman field h z from zero, we may cross several phase 
boundaries, e.g., from i-US-0 to US-2 to US-1, or from 
d-US-0 to US-1. Correspondingly, both kt and \zz are 
non-analytic at the phase boundaries (Fig. 2J: and0Jl), 
thus providing clear thermodynamic signatures of the 
topological quantum phase transitions described above. 

In summary, we investigated the effects of Zeeman 
fields and ERD SOC on the superfluidity of two- 
dimensional ultra-cold fcrmions from the BCS to the 
Bose regime. Wc constructed the ground-state phase di- 
agram of Zeeman field versus spin-orbit coupling, identi- 
fied bulk topological (Lifshitz) phase transitions between 
gapped and gapless superfluid. We also described the 
existence of a tri-critical point in the strong spin-orbit 
coupling regime. Lastly, we analyzed the excitation spec- 
trum and the momentum distribution, as well as thermo- 
dynamic quantities including the pressure, compressibil- 
ity, induced polarization, and spin-susceptibility, all of 
which provide signatures of the phase transitions. 

We thank Wei Zhang for discussions and ARO 
(W911NF-09-1-0220) for support. 
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